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We study a sample of approximately 16, 500 galaxies with Iacs,ab < 22.5 
in the central 38% of the COSMOS field (Scoville et al. (2006), this volume), 
which are extracted from a catalog constructed from the Cycle 12 ACS F814W 
COSMOS data set. Structural information on the galaxies is derived by fitting 
single Sersic (Sersic (1968)) models to their two-dimensional surface brightness 
distributions. 

In this paper we focus on the disk galaxy population (as classified by the Zurich 
Estimator of Structural Types, ZEST; Scarlata et al. (2006a), this volume), and 
investigate the evolution of the number density of disk galaxies larger than ap- 
proximately 5 kpc between redshift 2; ~ 1 and the present epoch. Specifically, 

we use the measurements of the half-light radii derived from the Sersic fits to 
construct, as a function of redshift, the size function $(r 1/2,-2) of both the total 
disk galaxy population and of disk galaxies split in four bins of bulge-to-disk 
ratio. In each redshift bin, the size function specifies the number of galaxies per 
unit comoving volume and per unit half-light radius ri/2- Furthermore, we use 
a selected sample of roughly 1800 SDSS galaxies to cahbrate our results with 
respect to the local universe. We find that: 

(i) The number density of disk galaxies with intermediate sizes (ri/2 ~ 5-7 kpc) 
remains nearly constant from ^ ~ 1 to today. Unless the growth and destruction 
of such systems exactly balanced in the last eight billion years, they must have 
neither grown nor been destroyed over this period. 

(ii) The number density of the largest disks (ri/2 > 7 kpc) decreases by a factor 
of about two out to 2; ~ 1. 

{iii) There is a constancy - or even slight increase - in the number density of 
large bulgeless disks out to 2; ~ 1; the deficit of large disks at early epochs seems 
to arise from a smaller number of bulged disks. 

Our results indicate that the bulk of the large disk galaxy population has com- 
pleted its growth by ^ ~ 1, and support that secular evolution processes produce 
- or at least add stellar mass to - the bulge components of disk galaxies. 

Subject headings: cosmology: observations — cosmology: large scale structure 
of universe — galaxies: formation — galaxies: evolution — galaxies: disks — 
galaxies: late- type — surveys 
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1. Introduction 

A detailed quantification of galaxy structure as a function of redshift is an essential 
tool to study the evolution of galaxies through cosmic time. In particular, the evolution of 
galaxy sizes is a direct diagnostic of their star formation and mass assembly history. For 
disk galaxies, several theoretical predictions have been made, which start being testable with 
observations: For example, under the assumption that the scale lengths of disks scale with 
their virial radii, their sizes are predicted to evolve as oc H{z)~'^/^ at a fixed halo mass 
(Fall & Efstathiou (1980), Mao et al. (1998)); Bouwens k Silk (2002), instead, assume a 
dependence of the star formation rate on the local gas density within the disk, and infall of 
metal-free gas, and predict the sizes of disks to evolve as R{z)/R{0) — 1 — 0.27 • z. 

Pioneering observations, however, generally appear to support a substantially weaker 
size evolution of the disk galaxy population than predicted by theory. Data from the Canada- 
France Redshift Survey (CFRS) for a few tens of large disk galaxies up to redshift 2; ~ 1 
indicate about a magnitude of surface brightness evolution for the disks, and a constant 
number density of large disks over this period (Schade et al. (1995); Schade et al. (1996); 
Lilly et al. (1998)). Similar results are reported by Roche et al. (1998). In contrast, Simard 
et al. (1999) interpret their data as indicating no surface brightness evolution in disks over 
the same redshift range. The work of Trujillo & Aguerri (2004) , on the other hand, provides 
further evidence for surface brightness evolution in galaxy disks, although a possible decrease 
in the sizes of disks by approximately 30% at 2; ~ 0.7 is also proposed as an alternative to 
explain the data. Using 1500 disk galaxies at redshifts 0.25 < z < 1.25 from the GOODS 
survey, Ravindranath et al. (2004) find at most a very modest decrease of the number of 
disks with half-light radius ri/2 > 4 kpc. The GEMS survey has been used to investigate the 
evolution with redshift of the luminosity-size relation for disk galaxies (Barden et al. (2005); 
see also Trujillo et al. (2005)), but no analysis of the evolution of the number density of disk 
galaxies as a function of galaxy size has been presented based on the GEMS data. Therefore, 
the observational evidence for a relatively constant number of (large) disks since 2; ~ 1 relies 
on very small samples extracted from small areas in the sky. 

In this article we use the COSMOS survey (Scoville et al. (2006), this volume) to inves- 
tigate how the number density of disk galaxies with sizes above 5 kpc - the threshold above 
which our sample is complete - evolves with cosmic time from 2; ~ 1 to 2; = 0. Thanks to 
the large statistics of COSMOS, we are able to extend the investigation so as to study this 

evolution as a function of the bulge-to-disk ratio of disk galaxies. Specifically, we use the 
ACS COSMOS database, and select a Iacsab < 22.5 sample from the 38% of the COSMOS 
area that was imaged by the HST during Cycle 12. The complete Iacs,ab < 22.5 sam- 
ple contains more than 16,000 galaxies. We perform two-dimensional GIM2D (Marleau & 
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Simard (1998)) fits to their surface brightness distributions, which we describe with a single 

Sersic law (Sersic (1968)). Free parameters in the fits arc: the total fiux Ffot, the Sersic 
index n, the half-light radius ri/2, the ellipticity e — 1 — b/a (with a and b the semimajor- 
and minor-axis), and the position angle 0. 

We compare the structural parameters obtained with the GIM2D surface brightness fits 
with the structural measurements derived using the Zurich Estimator of Structural Types 
(ZEST; Scar lata et al. (2006a), this volume). ZEST uses individual measurements of four 
non-parametric diagnostics of galaxy structure, namely, the concentration C, the asymme- 
try A, the Gini coefficient G, and the second order moment of the brightest 20% of the 
galaxy pixels, M20 (c.f. Scarlata et al. (2006a), this volume, and references therein), phis the 
elongation of the distribution of the light, to classify galaxies in three morphological types: 
early-types, disk galaxies or irregular and peculiar galaxies. The disk galaxies are further 
split by ZEST into four classes of bulge-to-disk ratio (B/D). At any redshift, the multi- 
dimensional ZEST classification grid substantially reduces the intermixture of structurally 
different galaxy populations compared with simpler approaches. The comparison of our 
parametric structural quantities with the diagnostics supplied by ZEST shows the expected 
correlations (e.g., between Sersic index n and concentration C) and, by doing so, illustrates 
the robustness of the results of our GIM2D fits. It also demonstrates that studies of selected 
galaxy populations which are classified by a simple cut in Sersic index n - as customarily 
done, for instance, in studies of the SDSS galaxy population and in other surveys of the 
high- 2; universe - are affected by mixing of galaxies with different structural properties. 

We use the ZEST classification to extract, from the total Iacs,ab < 22.5 sample, the 
subsample of approximately 12,000 disk galaxies (with different B/D ratios), and use the 
half-light radii provided by our Sersic fits to compute the number of disk galaxies per unit 
comoving volume and per unit half-light radius ri/2 as a function of redshift (the size function 
$(ri/2, z)). In order to compare the results at higher redshift with the local universe, we use 
a suitably selected sample of about 1800 SDSS galaxies, which have been "redshifted" to 
2; = 0.7 by properly re-scaling for fiux dimming, distance and image resolution (Kampczyk 
et al. (2006), in preparation). 

We find that the number density of disk galaxies with intermediate sizes ri/2 ~ 5-7 kpc 
remains nearly constant from 2; ~ 1 to today, while that of the largest disks (ri/2 > 7 kpc) 
shows a drop by a factor of approximately two. This deficit appears to stem mostly from a 
decline in the number density of bulged disks at early epochs.^ 



^In the rest of the paper, we consistently use AB magnitudes (Oke (1974)), and thus drop the subscript 
"AB" from our notation. We adopt Vim = 0.25, Qa + = 1 and a present day Hubble parameter of 
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2. Data analysis 

2.1. The initial catalog: 16,538 ACS-selected galecxies 

We use a sample of 16,538 COSMOS galaxies with I < 22.5, detected in the 259 ACS 
F814W images acquired during the HST Cycle 12 observing period (Scoville et al. (2006)). 
Details on the processing of the ACS data are given in Koekemoer et al. (2006); a discussion 
of the ACS-detected catalog used as the starting point to build our sample is presented in 
Leauthaud et al. (2006). In the following we give a brief overview over the steps leading to 
the production of the catalog: 

• The ACS-detected catalog was generated by running SExtractor (Bertin & Arnouts 
1996, version 2.4.3) twice on the reduced ACS images; the first time in a "cold" run 
with a configuration optimized for the detection of large, bright objects, and then in a 
"hot" run with a configuration optimized for small and faint sources. 

• The two resulting samples were then merged together to produce a final catalog by 
retaining all the "cold" detections plus the "hot" detections which fell outside the 
SExtractor segmentation map of any galaxy detected in the "cold" run. This ACS- 
based catalog contains approximately 55,700 galaxies down to / = 24. 

• Stars were removed by deleting the sources with SExtractor CLASS_STAR> 0.6 from 
the catalog. 

• The ACS-based catalog was inspected to remove all remaining over-deblended large 
galaxies as well as false detections from it. About 4.5% of the sources in the original 
ACS-based catalog were deleted from our final catalog after this visual check. 

The final and cleaned ACS-based catalog contains 55,651 galaxies down to / = 24, 
16,538 of which have / e [16, 22.5] and an ACS-stellarity parameter smaller than 0.6^. 
These roughly 16,500 galaxies constitute the sample that we study in this article. 



Ho = 70 km s"^ Mpc''^. 

^Notc that by applying the customary cut at a vahic of the stcUarity parameter of 0.9 some very bright 
stars would be included in the sample. We inspected the sources with stellarity larger than 0.6 to ensure 
that no galaxy was excluded from our sample. 
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2.2. Two-dimensional GIM2D fits 

The GIM2D (Galaxy IMage 2D) IRAF software package is designed for the quantitative 
structural analysis of distant galaxies (Marleau & Simard (1998), Simard et al. (2002)). 
GIM2D uses the Metropolis algorithm to converge to the analytical model, i.e. it carries 
out a Monte Carlo sampling of the likehhood function P{w\D, M), which measures the 
probability that the parameter combination w is the correct one given the data D and the 
model M. The Monte Carlo approach of samphng the complex multidimensional topology 
of parameter space has the virtue of converging consistently to the same best fit model for 
a wide range of initial guesses. 

In this analysis, we adopt a simple single Sersic (Sersic (1968)) profile to describe the 
two-dimensional surface brightness distribution of a galaxy's light: 



E(r) = E(ii!i/2) exp < ~k 



R 



1/2 



1/n 



(1) 



where the value of k is chosen such that R1/2 represents the radius containing half of the total 
flux. (Throughout this article we adopt the convention that R1/2 stands for the apparent 
half-light radius on the image, measured in pixels, while ri/2 denotes the physical half-light 
radius in units of kpc). The flexible form of the Sersic law has the advantage of parametrizing, 
through the variable exponent n, surface brightness distributions including the exponential 
radial fall-off of the light profile in bulgeless disks (n = 1), and the classical de Vaucouleur 
profile encountered in elliptical galaxies (n = 4) . 

In the case of a single Sersic model, GIM2D seeks the best fitting values for the following 
eight parameters: the total fiux Ffot (integrated to r = 00); the half-light radius R1/2', the 
ellipticity e = 1 — b/a, where a and b are the semimajor and semiminor axes of the brightness 
distribution; the position angle (p; the offsets dx and dy from the initially specified centre of 
the galaxy; the residual background level db and the value of the Sersic index n. 

The procedure followed to obtain Sersic models for all 16,538 galaxies in our sample 
involved the steps listed below: 

• Definition of the area to fit. For each galaxy, GIM2D was instructed (by the means of 
a mask image) to fit the distribution of light within an elliptical area with a semi- 
major axis equal to 1.5 Petrosian radn (Petrosian (1976))^ of the target galaxy and 



^The Petrosian radius rp is the radius at which the surface brightness in a thin annulus equals a given 
fraction of the average surface brightness within that radius, in our case a fraction of 20%. Since it is 
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an ellipticity equal to the flattening of the SExtractor segmentation map (the latter 
defines the pixels associated with the object). The orientation of this Petrosian ellipse 
is also provided by SExtractor. During the fitting, the algorithm considers both pixels 
assigned to the objects of interest and those flagged as background in the mask image. 
The goodness of the flt (i.e. the x^), on the other hand, is estimated based on the 
difference betweeen image and model for pixels within the Petrosian ellipse. 

• Generation of stamp images. For each galaxy wc extracted from the reduced ACS im- 
ages, a stamp image, centered on the coordinates of the ACS-bascd catalog. The stamp 
images were (i) no smaller than 10" x 10", and (ii) always sufficiently large to guar- 
antee a proper computation of the local background for all target galaxies. Point (ii) 
is ensured by choosing the stamp size such as to contain at least the same number of 
"background" pixels outside the masking lie within an ellipse with semimajor 
axis equal to a single Petrosian radius, centered on the target galaxy. 

• Generation of cleaned images. To optimize the quality of the flts, wc allowed G1M2D 
to improve the determination of the center of the galaxy with respect to the coordinates 
specified in the catalog. Furthermore, GIM2D requires and computes the value of the 
background to produce the fit to the galaxy light. Bright objects that are located 
near the galaxy to be fit can cause an inaccurate calculation of the background, and 
an incorrect estimate of the galaxy center. To prevent these potential errors from 
affecting the quality of our Sersic models, we cleaned the galaxy stamp images of 
all close companions. In this cleaning procedure pixels in the stamp image that were 
associated - based on a segmentation map accompanying the ACS-catalog - with sources 
other than the one of interest were replaced by values reproducing the properties of 
the surrounding background (c.f. Scarlata et al. (2006a), this volume, for details). We 
used these cleaned images in the next steps. 

• Symmetrization of galaxy images. Wc used G1M2D to produce a symmetrized version 
of the cleaned galaxy image. This was done by rotating the galaxy by 180°, subtracting 
the rotated image from the cleaned image, setting all pixels with values smaller than 
2a of the cleaned image's background in the difference image to zero, and, finally, 
subtracting this clipped difference image from the cleaned image. These symmetrized 
images were used as input to GIM2D in order to converge to a satisfactory model of 
the galaxy light distribution more reliably and quickly. 



computed from the ratio of two surface brightness values, the Petrosian radius is not as sensitive as other 
size measures to the effect noise has on objects at higher redshifts. It thus provides a robust and consistent 
definition of the surface over which to perform our GIM2D fits. 
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• Initial launch parameters and allowed parameter space. Given the robustness of GIM2D 
with respect to changes in the initial parameters (see Appendix A.l), we started all 
our fits with the same initial guesses of galaxy structural parameters, namely: / = 21, 
ri/2 = 10 px, e = 0.5, = 0° and n = 4. The calculation of a pure Sersic profile is 
achieved by fixing the ratio of fiux in the bulge to the total fiux, B/T, to unity. Struc- 
tural parameters were allowed to vary within the following ranges: / G [17.5,26.7], 
ri/2 e [1 pa;, 50 px], e G [0, 1], (f) G [-180°, 180°], n G [0.2,9]. We inspected the rela- 
tively small number of galaxies that appeared to genuinely require values outside these 
ranges, and allowed the parameters to converge to best fits outside those boundaries 
when our visual inspection confirmed the quality of the fits. 

• Adopted Point Spread Functions (PSFs). We select a different PSF for each individual 
object based on its position on the tile and the current focal length of the ACS due to 
in-orbit breathing at the time of the exposure (Rhodes (2006), this volume). 

2.3. Deletion of sources from the original ACS-based catalog 

A clear assessment of the selection effects affecting galaxy samples is essential to perform 
statistical studies of the evolution of galaxy populations. Sources had to be removed from 
our initial sample of 16,538 ACS-selected, Cycle 12 COSMOS galaxies with / G [16, 22.5] 
and ACS-stellarity smaller than 0.6 due to a variety of reasons: 

• A total of 465 sources (i.e., about 2.8% of the sample) were deleted from our final 
sample, due to the failure of GIM2D to converge. 

• 29 additional sources were also excluded on the basis of the of the G1M2D fits. 
After inspection of the G1M2D residuals and of the Xgim2D distribution, we fixed 
the threshold for deleting sources from our sample to Xgim2D — l^- While this may 
seem an exceedingly generous criterion, visual inspection of all GIM2D fits and their 
residuals for objects with Xgim2d > l-^ showed that, in the vast majority of the cases, 
GIM2D produces a fair model of the diffuse galaxy light, and the high vahics of Xgim2D 
in these objects were produced by structure such as knots of star formation, spiral arms 
or bars. Figure 1 and the first three entries in Table 1 display and list three random 
examples of objects with Xgim2D ^ [5, 15]. It is clear from the figure that the presence 
of spiral and other structure can lead to high values of Xgim2D even if the underlying 
symmetric component is reproduced well by the Sersic model. 

• Of the galaxies in our sample, 344 lacked ground-based photometric information in at 
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least five of eight bands; therefore, no accurate photometric redshifts could be derived 
for these sources. 

• A fourth selection criterion for the rejection of objects was apphed in part of the 
analysis, based on the value of xlhz supplied by ZEBRA (c.f. Section 4.1). As a test 
to evaluate the impact of high values of x^hz photometric redshift estimates, we 
adopted the approach of calculating the size functions (Section 4) both by applying no 

Xphz "^^t, and by deleting from mir sample sources with a Xphz above a threshold given 
by the value corresponding to two standard deviations in an ideal ^^-distribution of a 
varying number of degrees of freedom (which in our case are the number of filters for 
which ground-based photometric measurements are available). This cut - when applied 
- reduces the sample by 1369 objects. 

• Finally, objects with one ground-based set of coordinates but two distinct coordinates 
in the ACS-selected catalog must be expected to have ambiguous photometric infor- 
mation and they were therefore removed as well. A total of 271 objects fell in this 
category. 

3. The ZEST morphological classification of our COSMOS sample 

All COSMOS galaxies in our sample have been classified with ZEST, the Zurich Esti- 
mator of Structural Types, discussed in Scarlata et al. (2006a) in this volume. ZEST is a 
multi-dimensional classification grid that simultaneously uses: (i) the non-parametric diag- 
nostics of asymmetry A, concentration C, Gini coefficient G, and the second order moment 
of the brightest 20% of galaxy pixels M20, (e.g., Lotz et al. (2004)), plus (ii) the elongation 
of the galaxy's distribution of hght, to classify galaxy types. For each galaxy, the final ZEST 
classification provides quantitative information on galaxy structure, and a classification into 
either of the types 1, 2 or 3 for early-type, disk and irregular or peculiar galaxies, respectively. 
ZEST further splits the disk galaxies into four sub-classes of bulge-to-disk ratio, ranging from 
bulgeless disks (type 2.3) to bulge-dominated spirals (type 2.0). The Sersic indices presented 
in this paper were employed in ZEST to check and refine the four B/D subclasses of disk 
galaxies. 

Figures 2 to 7 show, for the different ZEST types and for disk galaxies with different 
B/D ratios, a few examples of GIM2D fits picked randomly among representatives of the 
morphological class in question. In each figure we plot (from left to right) the original 
stamp image, the cleaned image, the GIM2D model and the fractional residuals between the 
observed galaxy and the analytical model. It is clear from these examples that early-type 
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and disk galaxies are well described by the analytical fits with the Sersic profile. They do 
well in reproducing the smooth component underlying the galactic light and, as is visible in 
the residuals, leave behind only the sub-structure within the galaxies, such as spiral arms 
and bars. On the other hand, the GIM2D fits for the irregular galaxies are understandably 
less reliable, as the smooth and symmetric analytical models cannot capture the richness of 
structure of this class of objects. 

To further highlight the robustness of the GIM2D fits, we show in Figure 8 the com- 
parison of our values of the Sersic index with the four non-parametric diagnostics of galaxy 
structure computed to derive the ZEST classification of the galaxies in our sample. The 
expected correlations between diagnostics such as concentration and Sersic index, as well as 
the tightness of these relationships, are an independent confirmation of the reliability of the 
analytical fits. 

The left hand side of Figure 9 shows the distribution of Sersic indices n as a function 
of ZEST galaxy type. The figure illustrates that, as expected, the bulk of late-type galaxies 
(including the type 3 irregular galaxies) typically shows low values of n and preferentially 
lies close to the typical value n = 1 of disk-like exponential profiles. The distribution of 
early-type galaxies, on the other hand, reaches its peak at n ~ 3.5. The right panel of the 
figure shows the distribution of n for disk galaxies of different B/D ratio; by construction, 
the mode shifts to ever higher values of n as the bulge component gains in importance. It is 
evident from Figure 9 that a separation of galaxy types according to a simple cut in Sersic 
index brings about a significant mixture of different galaxy populations. Depending on the 
scientific goals, however, it may be important to use samples as uncontaminated as possible 
by other galaxy types. In the following we thus adopt the more precise ZEST classification 
of galaxy types - rather than a cut in Sersic index - to investigate the evolution of the size 
function of disk galaxies (Section 4.3). 

As a final illustration of the COSMOS sample under investigation, we show the distribu- 
tions of physical sizes ri/2 in kpc in Figure 10. ZEST types and sub-types are displayed with 
different line styles, as in Figure 9. The distribution of half-fight radii of disk-dominated 
galaxies (morphological type 2) peaks at roughly 2.75 kpc, that of eUiptical galaxies (type 1) 
around 2 kpc. In the right panel of the figure, the peaks of the distribution of half-light radii 
shift to progressively smaller sizes with increasing B/ D-ratio, as is expected for objects with 
successively more centrally concentrated surface brightness profiles. The only exception to 
this trend are the most bulge-dominated sources among the disk galaxy population (type 
2.0), which break out of this sequence. The ZEST class 2.0 is defined as containing objects 
with a spheroidal component (and Sersic index) similar to the type 1 early-type galaxies, 
but, as opposed to the latter, showing an obvious disk component as well. For many scien- 
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tific applications, it might prove beneficial to merge this morphological class with the type 1 
early-type galaxies (c.f., for instance, Scarlata et al. (2006b) for a discussion on the evolution 
of early-type galaxies in COSMOS, in which the effects of including the ZEST type 2.0 class 
in the sample is discussed in detail). Here we emphasize that the / < 22.5 COSMOS sample 
under scrutiny encompasses, for each galaxy type individually, a broad range in physical 
sizes and allows for statistically sound statements thanks to the large number of available 
objects. Thus, with proper care given to the effects of incompleteness, this permits us to 
study the evolution with cosmic time of the comoving number density of galaxies of a given 
size. 

In this analysis we focus exclusively on the disk galaxies - classified as such by ZEST 
- and postpone the investigation of other galaxy populations to future publications. Of 
the 14,520 galaxies reliably fitted with GIM2D that remain in the sample after applying the 
selection criteria in magnitude, stellarity parameter, Xgim2Dj X^hz photometric reliability 
(see Section 2.3), 11,744 galaxies are classified as disk galaxies on the ZEST COSMOS 
classification grid. If no cut based on the value of Xphz applied, the remaining sample 
of disk galaxies contains 12,701 objects. Figure 11 shows the distribution with redshift of 
the galaxies in our sample and Tables 2 and 3 summarize in more detail how the disks and 
other galaxy types are assigned to the different investigated redshift bins up to z ~ 1. This 
information is provided both for the case when a cut based on the value of Xphz applied 
to the sample, and the case when no such cut is apphed. We have tested that the results 
obtained with and without the cut in Xphz show no significant difference. We therefore show 
our results using the sample with the cut in Xphz- Doing so implies the exclusion of a total 
of 8% of the original ACS-based 7 < 22.5 Cycle 12 COSMOS sample. 



4. The Evolution of the Size Function of disk galaixies in COSMOS 

4.1. ZEBRA photometric redshifts 

Several codes and approaches have been used to derive accurate photometric redshifts for 
COSMOS galaxies down to faint magnitudes (c.f. Mobasher et al. (2006) in this volume for 
a review). In our work, we adopt the photometric redshift estimates obtained by Feldmann 
et al. (2006) for the COSMOS galaxies of our ACS-selected / < 22.5 galaxy sample. These 
are computed with the Zurich Extragalactic Bayesian Redshift Analyzer (ZEBRA). 

ZEBRA produces two separate estimates of the photometric redshifts of individual 
galaxies: A maximum likelihood estimate and a fully two-dimensional Bayesian estimate 
in the space of redshift and templates. In both approaches, ZEBRA uses an iterative tech- 
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nique to automatically correct the original set of galaxy templates to best represent the 
spectral energy distributions of the galaxies in different redshift bins. The availability of a 
"training set" of spectroscopically derived zCOSMOS redshifts (c.f. Lilly et al. (2006), this 
volume) for a small fraction of the whole considered photometric sample allows for a precise 
calibration of the ZEBRA photometric redshifts and thus for an optimal correction of the 
galaxy templates. 

The ZEBRA uses the available ancillary ground-based multi-band photometry to esti- 
mate the photometric redshifts of the COSMOS galaxies. B, V, g, r', i', and z' fluxes are 
measured with Subaru (Taniguchi ct al. (2006)), u* fluxes with CFHT and the 4 m telescopes 
at Kitt-Pcak and CTIO provide the Ks photometry (Capak et al. (2006)). 

In the following, we use the ZEBRA maximum likelihood photometric redshifts of Feld- 
mann et al. (2006) for our discussion of the evolution of the size function of disk galaxies. 
The maximum likelihood ZEBRA photometric redshifts have, at all redshifts of interest, an 
accuracy of Az/{1 + z) ~ 0.03 if compared with the zCOSMOS spectroscopic redshifts (Lilly 
et al. (2006)). The application of ZEBRA to our sample shows some dependence of the pho- 
tometric redshifts on whether small systematic offsets (of order 0.05 magnitudes or smaller), 
which are detected by the code, are apphed to the calibration of the Subaru data. This has, 
however, no substantial impact on our analysis, as we show in the following by using, when 
relevant, the two ZEBRA estimates of the photometric redshifts obtained with and without 
corrections for these photometric offsets. The inclusion or exclusion of bad ZEBRA fits from 
our final sample does not change the conclusions of this manuscript. 



4.2. Derivation of the Size Function and assessment of completeness 

In each redshift bin, the size function $(ri/2,-z) measures the number of galaxies per 
unit comoving volume and per unit half-light radius ri/2- In order to chart the evolution of 
the number density of disk galaxies with a given size as a function of redshift up to 2; ~ 1, we 
split our sample of disk galaxies into four different redshift bins of width Az = 0.2, centred 
on the redshifts 0.3, 0.5, 0.7 and 0.9. 

There are two different considerations in constructing the size function which should be 
clearly distinguished. The first is the variation of sampling volume for different individual 
objects that are included in the sample. This can be fully taken into account using the 
Knaa:-formalism (Schmidt (1968); Felten (1976)) in which each galaxy is weighted with the 
reciprocal of the comoving volume it could occupy while still satisfying the selection criteria 
of the sample. For bright objects this volume will be larger, thus leading to a smaller 
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weight and correcting for their relative rarity with respect to fainter members of the galaxy 
population: 

$(ri/2,^)xAr,^/^ = E^-^. (2) 

The summation is carried out over all sources 7 with a photometric redshift within the redhift 
bin in consideration. Using that the comoving distance D(z) in a ACDM cosmological model 
with fl\ — 0.75 and fi^ = 0.25 is given by 

the volume Vj^ax is calculated for each source according to: 

dz 



Vm.ax ^ ^ 



(4) 



Here dO, is the effective sohd angle covered by the survey during Cycle 12 which is 2734 



arcmin^. 



The lower and upper bounds of integration, Zmin and Zmax, define the redshift range 
within which any particular galaxy in the sample would have entered the sample. It is 
constrained in three ways: (a) the minimum and maximum redshifts of the relevant redshift 
bin; (b) the minimum and maximum redshifts at which the sample magnitude limits are 
satisfied, and (c) the minimum and maximum redshifts at which the size limits in the sample 
are fulfilled. The size limit of relevance is that the object must have sufficient pixels to be 
classified by ZEST (FWHM > 0.15 arcsec) and can be calculated from the condition that 
the galaxy have the same physical size ri/2 oc Da ■ Ri/2[px] at all redshifts, regardless of 
its observed scale (RAi^iJipi-siiA — —51/2 — y practice the size constraint is the least 
significant because of the small variation in angular diameter with redshift at high redshifts 
and because of the high angular resolution of the HST images. Of more importance are the 
limits in redshift imposed by the magnitude cuts. In the case of a bounding magnitude of 
Iiim - in the present analysis either 7 = 16 or 7 = 22.5 - the requirement of identical absolute 
magnitude leads to the condition: Mb — him — 5 log ^^^1^^^ + Kfsiaw,b{z)^ which must 
be solved for the limiting redshift Zum'^. 

There may also be an observational surface brightness selection. It is, however, unim- 
portant in practice when compared with the surface brightness constraint induced by a 



^Here Kp&uw, b denotes the K-correction as defined through the equation Iab = Mb + 
5 log {Dl{z)/IQ pc) + Kj^b (cf. Hogg et al. (2002)). 
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selection in absolute magnitude. The difficulty is that the size function is only partially 
sampled, because of the magnitude limits of the sample. Thus the size function constructed 
above represents the integral of the bivariate size-luminosity function above some limiting 
luminosity. We adopt the approach of Lilly et al. (1998) to deal with this problem. We (a) 
limit the sample in equation (2) to objects lying within a specified range of luminosity, and 
allow this luminosity range to change with redshift according to the expected luminosity 
evolution of individual galaxies. The integration over luminosity of the bivariate function is 
thus explicit, albeit model dependent. Secondly, we consider, (b), at each redshift, the limit- 
ing size above which one can be fairly confident that the sample is more or less "complete" , 
i.e. that there are rather few objects of this size, or greater, lying below the magnitude limit 
because of their low surface brightness. Of course, luminous low surface brightness galaxies 
are known to exist and the sample can never be absolutely complete, but we may assume 
that the size function above this limiting size will be a good approximation to reality and it 
will, in any case, represent a lower limit to the actual number density. It is important to ap- 
preciate that this apparent surface brightness selection arises from the inevitable selection in 
absolute magnitude (at a given size) rather than any direct observational surface brightness 
selection (c.f. Lilly et al. (1998) for a detailed discussion of this matter). 

Figure 12 shows the distribution of Mb and ri/2 in the sample as a function of redshift. 
As discussed in the previous paragraph, we limit the set of galaxies to those above a redshift- 
dependent magnitude limit of Mb — (—19.6 — 2;). It refiects an assumed luminosity evolution 
by about one magnitude due to a passively aging stellar population since 2; ~ 1. The 
horizontal dotted lines in Figure 12 illustrate the effect of the cut as it would be applied 
at the central redshift of the different bins. Their loci change from an absolute magnitude 
of Mb = —20.5 at redshift 0.9 to Mb = —19.6 at redshift zero. In Figure 12 the inclined 
dotted lines are lines of constant surface brightness, shifted according to the assumed surface 
brightness evolution. The figure illustrates that our sample is fairly complete above ri/2 ~ 5 
kpc, i.e. the bulk of galaxies with sizes above this threshold lies above the horizontal line at 
all redshifts. We therefore limit our analysis to galaxies with half-light radii exceeding this 
size. 



4.3. Errors 

Sizes, magnitudes and photometric redshifts are all subject to uncertainty. In order 
to assess the resulting inaccuracy of the size function $(ri/2,2;), we varied these quanti- 
ties according to a Gaussian distribution centred on the measured values and having an 
appropriate standard deviation. In the case of the photometric redshifts, for instance, the 
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standard deviation rises with increasing redshift as ~ 0.03- (1 + 2;). For the magnitude, we 
used the standard deviation o"/ given in the ACS-COSMOS catalog. Errors for the sizes are 
tabulated for each galaxy by G1M2D in the form of upper and lower 99% confidence levels, 
corresponding to ScTi^^^^ for the half-light radii. 

Prom the results of 101 test runs with errors introduced on all three of the above 
quantities, we calculated the resulting scatter for each bin of the size function. The size 
of the errors given on the grounds of these tests corresponds to the interquartile range of 
the 101 values at every measurement point on the curve $(ri/2,-2)- The interquartile range 
spans 50% of a data set, and eliminates the influence of outliers by effectively removing the 
highest and lowest quarters of the data range. The distance from the median of the 101 error 
realisations to the first and third quartile give the lower and upper error bar respectively. 
In the plots of the size function (c.f. Figures 13 and 14) these asymmetric errors have been 
reported in each individual data point and the area between them shaded in grey. 

The Poisson noise in each bin of the size function (owing to the statistical manner in 
which each bin is populated) remains, however, the dominant source of random error in the 
size function. The Poisson error is given by: 



with the summation carried out over all objects i in a given bin. Poisson errors are reported 
symmetrically around the data points in each bin of the size function. 

Systematic uncertainties in the GIM2D measurements of the half-light radii (discussed 
in Appendix Al) could affect the slope of the size functions presented in Figures 13 and 14 
for the COSMOS galaxies, and in Figures 16 and 17 for our comparison sample of SDSS 
galaxies (section 4.4 below). It should be noted that any systematic errors arising from 
the fits should be identical at 2; = and z — 0.7 and thus the SDSS-normalized densities 
of Figures 18 and 19 - on which we base our discussion and conclusions - are free of such 
systematics. 



To quantitatively assess the redshift evolution of the number density of disk galaxies of a 
given size (and as a function of B / D-ratio) , it is important to realize that a large overdensity 
in the COSMOS field at z < 0.4 affects the number counts in the lowest COSMOS redshift 
bin. Therefore, to carry out our comparison and determine how the number of disk galaxies 




(5) 



4.4. Comparison with the local SDSS galcixy population 
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with scalelengths larger than about 5 kpc changes with redshift from 2; ~ 1 to z ~ 0, we 
repeated our analysis on a set of 1876 galaxies extracted from the SDSS Data Release 4. 
To fully include all observational effects, the images of the SDSS galaxies were degraded 
so as to appear as they would in the COSMOS images at z = 0.7. The generation of 
these simulated images will be described in detail elsewhere (Kampczyk et al. (2006)). To 
summarize, the SDSS galaxies were selected from the redshift range z G [0.015, 0.025]. Their 
g-hand images were then transformed to how such galaxies would appear in the F814W 
COSMOS ACS images if they were to he a,t z — 0.7, at which redshift the passbands are well 
matched. Thus, this "redshifting" of the SDSS galaxies needed only to take into account 
the different pixel scales, and point spread function, and the cosmological surface brightness 
dimming. Two sets of images were produced, one without surface brightness evolution and 
one with an assumed brightening to high redshift of Afi = z magnitudes arcsec"^. (We will 
henceforth refer to these two samples as "unbrightened" and "brightened", respectively). No 
size evolution was considered. The simulated galaxies at z = 0.7 were then added into the 
COSMOS ACS images to reproduce also the same issues of image crowding and noise. These 
z — 0.7- "simulated" SDSS galaxies were then analyzed following the identical procedure 
adopted for the real COSMOS objects. Stamps sized 10" x 10" or more were extracted for 
all galaxies. As before, the original images were cleaned if necessary and then symmetrised 
prior to performing the GIM2D fits. After applying exactly the same selection criteria with 
respect to magnitude and Xgim2Dj 457 and 1294 objects remained in the unbrightened and 
brightened samples, respectively. The distribution of absolute B-band magnitudes, half- 
light radii and Sersic indices among the two sets of artificially redshifted SDSS galaxies is 
compared to that of COSMOS galaxies in the range z G [0.6, 0.8 [ in Figure 15. 

Absolute B-band magnitudes for the SDSS galaxies were derived by determining the best 
fit spectral energy distribution template with ZEBRA, and then calculating the i^'-correction 
for that spectral type, and the corresponding distance modulus, at redshift z ~ 0.7. The 
SDSS sample of Kampczyk et al. (2006) was extracted as a volume limited sample from the 
local Universe; its equivalent volume was obtained by computing the luminosity function for 
the 1876 SDSS galaxies, and setting it equal to the global SDSS luminosity function. We 
used the so-derived equivalent volume to compute the size function for the SDSS sample. 



5. Results and discussion 

Figures 13 and 14 show the size functions for the total disk galaxy population and 
for disk galaxies of different S/D-ratios, derived as a function of redshift above an evolving 
absolute magnitude hmit Mb — (—19.6 — 2;) (c.f. Section 4.2). As discussed before, these size 
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functions should be most complete for sizes larger than roughly 5 kpc. Above this threshold 
we have thus fit them with straight lines of log{^) vs. ri/2, i.e. with an exponential size 
function. The fits of the aforementioned functional form were calculated for two sets of 
photometric ZEBRA redshifts; one obtained without (solid line) and the other with (dashed 
line) corrections to the photometric catalog (with no significant change to our conclusions). 
The reported size function is based on the uncorrected photometry and the parameters of 
the corresponding solid black lines are given in Table 4. The size function of disk galaxies 
as a whole appears remarkably constant with redshift across the range 0.2 < z < 1.0. We 
detect, however, a weak trend for the slope of the size function to steepen with redshift, a 
result which applies to the global disk population and also separately to all sub-classes of 
bulge-to-disk ratio. 

These findings are strengthened by comparing directly the COSMOS size function in the 
range z e [0.6, 0.8[ with that derived from the local SDSS sample ( "redshifted" to z ~ 0.7). 
In Figures 16 and 17, the SDSS-based size functions are plotted together with the COSMOS 
size functions in the bin z G [0.6,0.8[. The size function of the brightened SDSS sample 
is reported with black points and error bars (and labeled as 'SDSS*'); in grey we show the 
unbrightened one ('SDSS'). The difference is due to the brightening of the galaxies which, 
in combination with the luminosity cut at Mb ~ —20.3 at z = 0.7 (c.f. Figure 12) leads to a 
greater number of smaller objects in the brightened set of SDSS galaxies. Given the evidence 
for luminosity evolution in all classes of galaxies, we would expect the brightened sample to 
provide the best comparison and in Figure 15 the latter indeed matches the distributions of 
Mb, n and ri/2 of COSMOS galaxies best. The line parameters of a hnear regression on the 
SDSS size function in the same range of sizes as in the COSMOS data set are listed in Table 
5. Figures 16 and 17 show a striking similarity in the shape of the size function over half 
the Hubble time. Furthermore, the steepening of the size function in COSMOS relative to 
locally is confirmed in that there is a small deficit of the largest disks (ri/2 > 7 kpc) in the 
COSMOS sample compared to the artificially redshifted local sample. 

To further investigate this issue we show in the two following figures the density of disk 
galaxies in COSMOS normalized to the local size function of (brightened) SDSS sources for 
two size bins, for all disks (Figure 18) and differentiated by bulge to disk ratio (Figure 19). 
In Table 6 we list the relative densities p = ^cosmos/^sdss* for disk galaxies as a function 
of redshift and B / D-c&Xegoij. Here ^cosmos denotes the comoving number density per unit 
size, as presented in Figures 13 and 14 for the COSMOS sample in the studied redshift bins. 
^SDSS* stands for the (artificially redshifted and brightened) SDSS density in the 0.6 — 0.8 
redshift bin (as presented in Figure 16 for the entire disk galaxy sample, and in Figure 17 for 
the sample split in the four B/D ratio bins). In the table, p< stands for the relative density 
of disks with intermediate sizes of ri/2 £ [5 kpc, 7 kpc[ and p> for that of large disks with 
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f\/2 £ [7 kpc, 17 kpc[. At 2; = 0, the normalized density point ^sdss*/^sdss* — 1 has been 
plotted as well, so as to show the error bar on the SDSS z = density estimate. Within 
the uncertainties (and excluding the lowest redshift bin from COSMOS because of the large 
cosmic variance), there is no evidence for a variation with redshift in the size function of 
disks in the range ri/2 G [5 kpc, 7 kpc[. Disk galaxies of these intermediate sizes have neither 
grown nor been destroyed from 2; = 1 to 2; = or, alternatively, the growth and destruction 
of such systems exactly balanced over this period. The larger disks, however, show a weak 
decline by a factor of nearly two over the same time span, even though the sample should 
be more complete for these larger disks as discussed before. 

Figures 14, 17 and especially Figure 19 and the values of the slopes a given in Table 4, 
consistently show the steepening of the slope of the size function to be the more pronounced 
for galaxies with higher B / D-raiios. This suggests that the deficit of large disks at ear her 
epochs is primarily due to a deficit of bulged disk galaxies. In principle disks with the highest 
bulge-to-disk ratio could be classified as early-types by ZEST and thus excluded from our 
sample, yet the direct comparison with the SDSS sample - which was classified in precisely 
the same fashion - argues against this interpretation. Moreover, the stated result remains 
valid even if type 2.0 galaxies are omitted from the discussion. Morphological /^-corrections 
could play a role in the observed disappearance of about 60% of the large disks with bulges, 
due to color differences between bulges and disks; we are studying this issue in more detail. 

The constancy - or even slight increase - of the number density of hulgeless disks out 
to 2; ~ 1 is at variance with models based on hierarchical formation or infall scenarios 
(e.g., Mao et al. (1998); Bouwens & Silk (2002)). On the other hand, a constant or slightly 
increasing number density of large bulgeless disks at 2; ~ 1, in combination with a decrease 
of similarly sized bulged disks to this redshift, supports a scenario where secular evolution 
processes form or at least grow the sizes of bulges at late stages in the history of the universe. 
Independent evidence for bulge-building secular evolution within disk galaxies comes from 
numerical experiments (e.g., Norman et al. (1996); Debattista et al. (2004, 2005, 2006)), 
observations of nearby "young" bulges with disk-like stellar density profiles (see Kormendy 
& Kennicutt (2004) and Carollo (2004) for reviews and references; Carollo (1999), Carollo 
et al. (1998, 2001, 2002), Carollo (2006)), and also from the lack of evolution in the fraction 
of large-scale bars in disk galaxies since z — 1 (Sheth et al. (2003); Jogee et al (2004)). 
Our work provides strong evidence that, by 2; ~ 1, the disks have not only grown to their 
current sizes, but also that they have become massive and stable enough to allow for large- 
scale internal dynamical instabilities to take place and to grow bulges inside disks, thereby 
shaping the Hubble sequence that we observe today. 
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6. Conclusions 

We have used measurements of the half-hght radii obtained with analytic fits to the 
surface brightness distribution of roughly 12,000 disk galaxies in COSMOS to chart the 
evolution of the number density of intermediate-sized and large disks since redshift z ~ 1. 

Our analysis has shown a general constancy of the size function for disk galaxies with 
half-light radii larger than about 5 kpc. This constancy, initially proposed by (Lilly et al. 
(1998)) based on much smaller statistics, and here clearly demonstrated from a much larger 
sample of disk galaxies of any bulge-to-disk (B/D) ratio, suggests that the massive disk 
galaxy population was largely in place with its current properties when the universe was 
only half as old as today. 

Among the observed general constancy of the size function there is evidence for some 
subtler and possibly quite important trends: (i) The number density of disk galaxies of all 
B/D-ratios in the intermediate size range ri/2 £ [5 kpc, 7 kpc[ remains constant to a good 
approximation up to 2; ~ 1. (ii) The number density of the largest disks with half-light radii 

£ [7 kpc, 17 kpc[ at 2; ~ 1 appears however to be only 60% of today's value in the disk 
galaxy population as a whole, (iii) When this population of large-sized disk galaxies is split 
by the values of their B/D-ratio we find that the deficit at high redshift stems from disk 
galaxies with a bulge component. It is conceivable that color differences between bulges and 
disks contribute to this result and we will study their effect in more detail. Bearing this in 
mind, however, it is nevertheless likely that we might be witnessing a secular transformation 
of bulgeless disks into disks with a bulge component. 
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A. Reliability of the GIM2D fits 

We have conducted several tests to assess the reUabiUty of the GIM2D fits. It must be 
emphasized that even objects with a high value of Xgim2D represent successful fits to the 
underlying symmetric component of the surface brightness profile in most cases. The 
value supphed by GIM2D is the sum over all pixels within the fitted area of the squared 
difference between model and observed data. Assuming that the pixel noise is purely due 
to Poisson photon statistics, this sum is then normalised to the variance of the pixel noise. 
Figure 1 clearly shows that high Xgim2D values mostly reflect the presence of spiral arms 
and other features of galactic structure that cannot be captured by the single-Sersic models, 
which continue however to provide a good description of the underlying galaxy light. 

A.l. Simulations 

To quantitatively test the robustness of the GIM2D flts, we performed a set of sim- 
ulations of an equal number of pure exponential disks (n = 1) and pure de Vaucouleur 
galaxies (n = 4) with a range of half-light radii and magnitudes (and ellipticities) covering 
the observed parameter space. 

It should be emphasized that the vast majority of disk galaxies in our sample is well 
described by em — 1 Sersic profile (c.f. Figure 9) and that none of the (disk-selected) galaxies 
in our sample are well represented by a smooth pure de Vaucoleur profile which would be 
appropriate for a disk-less spheroid. By construction, the ZEST-selected disk galaxies of our 
study posess a visible disk component, even if the single-Sersic fits that we have performed 
return an n = 4 profile. Nonetheless, we include the extreme pure de Vaucouleur cases in our 
simulation sample, to derive conservative estimates for the systematic errors in our GIM2D 
measurements. 

The simulated galaxies were generated with the IRAF task mkobjects, which automat- 
ically adds a component of Poisson noise to the flux of the source. Source-independent 
background noise was estimated from flve regions of pure sky in different ACS-tiles of the 
COSMOS survey. An average standard deviation from these regions was used to scale a 
Gaussian white noise distribution, in order to match the average background characteristics 
of the COSMOS images. The Sextractor measurements and the GIM2D fits for the simulated 
galaxies were carried out in an analogous way to the measurements and fits of the COSMOS 
galaxies (c.f. Section 2). 

Figures 20 and 21 summarize the main results derived from the simulated galaxies. Fig- 
ure 20 displays the distribution of normalized difference between the input and the recovered 
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Sersic indices n, as a function of galaxy magnitude and separately for the n — 1 and the 

n = 4 galaxies. Figure 21 shows the comparison between the input half-light radii and those 
derived from the G1M2D fits; in particular, the figure shows the distribution of the normal- 
ized difference between G1M2D and input half-light radii - separately for the n = 1 and the 
n = 4 galaxies as in the previous figure - as a function of magnitude (upper half) and of 
galaxy size (lower half). 

The GIM2D fits tend to underestimate the galaxy sizes. There is a dependence on 
galaxy magnitude and size, and on the shape of the profile: For example, the radii are 
underestimated by 3% and 10% for n = 1, Iab > 21 galaxies with half-light radii smaller or 
equal and larger than 10 pixels, respectively. The similar comparison for the n = 4 galaxies 
yields systematic uncertainties of order 20% and 40%, respectively. 

These systematic effects could affect the slopes of the size functions presented in Figures 
13 and 14 for the COSMOS galaxies, and in Figures 16 and 17 for the comparison sample 
of SDSS galaxies. We are, however, confident that these systematic errors on the GIM2D 
half-hght radii should not affect our conclusions. First, an influence of such systematics on 
the measured size function would actually cause an increase of the number density of large 
bulgeless systems at z ~ 1. Second, and most important, the approach of using a sample of 
SDSS galaxies redshiftcd to z=(}.7 to anchor the COSMOS size function to the local values 
ensures that these potential systematic errors on the size measurements for any disk galaxy 
type - including the most centrally concentrated, largest and/or faintest disk galaxies in our 
sample - should be identical at 2; = and z = 0.7. 



A. 2. Susceptibility to PSF Variations 

Although we did not expect the outcome of single Sersic fits to depend strongly on small 
variations of the PSF, we carried out tests with varying PSFs on a few COSMOS galaxies. 
The outcome of such a test for one of these objects is illustrated in Figure 22. The six 
panels show the distribution of parameters put out by G1M2D if the galaxy is fitted with 
250 different PSFs extracted in a (spatially) regular fashion from the total of 2500 PSFs 
generated by Rhodes (2006) for each ACS tile to study weak lensing maps. The results of 
the PSF tests confirm that the GIM2D fits are very robust to possible variations of the PSF 
as occur due to spatial variations across the ACS-tiles, for changes in the Sersic index and 
half-hght radius are well below the 10% level. 

We thank the anonymous referee for his valuable comments and suggested improvements 
to the original manuscript. 
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Table 1. Physical properties of the galaxies shown in Figures 1 to 7". 
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" The galaxies are listed in the order in which they appear in Figures 1 to 7, 
i.e. the first galaxy listed is the uppermost one in Figure 1, the third one the 
lowermost object in Figure 1, the fourth entry the uppermost galaxy in Figure 2, 
and so forth. 
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Table 2. Number of galaxies as a function of redshift when a 2cr-cut in Xphz applied 



rcdshift total sample type 2.3 type 2.2 type 2.1 typc2.0 type 1 type 3 
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Table 3. Number of galaxies as a function of redshift without a 2(7-cut in x^hz 



redshift total sample type 2.3 type 2.2 type 2.1 type2.0 type 1 type 3 



0.2 < z < 0.4 
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248 


0.8 < z < 1 


2840 


800 


552 


415 


287 


507 


270 


1 < z < 1.2 


661 


188 


134 


96 


52 


93 


96 
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Table 4. Line parameters of the COSMOS size function for radii ri/2 > 5 kpc". 



redshift 


type 


a 


/? 




0.3 


type 2 


-0.200 ± 0.021 


-2.278 ± 0.167 


0.912 


0.5 


type 2 


-0.215 ± 0.024 


-2.494 ± 0.180 


0.834 


0.7 


type 2 


-0.240 ± 0.017 


-2.225 ± 0.125 


5.363 


0.9 


type 2 


-0.294 ± 0.031 


-1.891 ± 0.215 


0.311 


0.3 


type 2.3 


-0.253 ± 0.070 


-2.602 ± 0.501 


0.204 


0.5 


type 2.3 


-0.221 ± 0.071 


-2.910 ± 0.485 


0.279 


0.7 


type 2.3 


-0.283 ± 0.034 


-2.321 ± 0.237 


2.940 


0.9 


type 2.3 


-0.291 ± 0.049 


-2.130 ± 0.331 


0.800 


0.3 


type 2.2 


-0.234 ± 0.039 


-2.372 ± 0.293 


0.457 


0.5 


type 2.2 


-0.236 ± 0.039 


-2.735 ± 0.291 


1.641 


0.7 


type 2.2 


-0.204 ± 0.022 


-2.913 ± 0.171 


6.963 


0.9 


type 2.2 


-0.260 ± 0.035 


-2.696 ± 0.254 


4.164 


0.3 


type 2.1 


-0.149 ± 0.061 


-3.468 ± 0.500 


0.088 


0.5 


type 2.1 


-0.171 ± 0.050 


-3.710 ± 0.446 


0.060 


0.7 


type 2.1 


-0.198 ± 0.043 


-3.446 ± 0.330 


1.562 


0.9 


type 2.1 


-0.214 ± 0.046 


-3.528 ± 0.377 


0.155 


0.3 


type 2.0 


-0.125 ± 0.029 


-3.585 ± 0.303 


0.981 


0.5 


type 2.0 


-0.129 ± 0.031 


-4.031 ± 0.326 


2.274 


0.7 


type 2.0 


-0.218 ± 0.046 


-3.304 ± 0.354 


0.082 


0.9 


type 2.0 


-0.242 ± 0.068 


-3.458 ± 0.533 


0.136 



Line parameters and of the linear regressions of the form 
log(^) = a ■ + /3 applied to the size functions of Figures 13 and 
14 at sizes ri/2 >5 kpc (solid black lines in the figures). The values 
of slope a and y-axis intercept /3 are given for the morphological 
class type 2 as a whole, and for the diflferent B/D-ratio sub-classes. 
See Figures 13 and 14 for the graphical representation of the size 
function derived from the COSMOS galaxies. 
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Table 5. Line parameters of the size function derived from artificially redshifted SDSS 

galaxies for radii ri/2 > 5 kpc". 



sample 


type 


OL 


/3 




SDSS 


type 2 


-0.159 ± 0.027 


-3.051 ± 0.234 


2.589 


SDSS 


type 2.3 


-0.172 ± 0.056 


-3.449 ± 0.494 


0.429 


SDSS 


type 2.2 


-0.110 ± 0.034 


-3.877 ± 0.341 


1.942 


SDSS 


type 2.1 


-0.128 ± 0.080 


-4.013 ± 0.675 


1.351 


SDSS 


type 2.0 


-0.100 ± 0.304 


-4.158 ± 2.046 


0.002 


SDSS* 


type 2 


-0.197 ± 0.026 


-2.470 ± 0.201 


0.939 


SDSS* 


type 2.3 


-0.176 ± 0.040 


-3.060 ± 0.326 


2.434 


SDSS* 


type 2.2 


-0.217 ± 0.055 


-2.727 ± 0.404 


0.402 


SDSS* 


type 2.1 


-0.107 d= 0.050 


-4.031 ± 0.503 


0.244 


SDSS* 


type 2.0 


-0.142 ± 0.048 


-3.743 ± 0.431 


1.683 



^ Line parameters and of the linear regressions of the form 
log{^) = a ■ + P applied to the size functions of the bright- 
ened (denoted as 'SDDS*' in the table) and unbrightened (denoted 

'SDDS') artificially redshifted SDSS galaxies at sizes ri/2 >5 kpc. 
The values of slope a and y-axis intercept /3 are given for the mor- 
phological class type 2 as a whole, and for the different B/D-ratio 
sub-classes. See Figures 16 and 17 for the graphical representation 
of the size function derived from the SDSS galaxies. 
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Table 6. Relative densities p — ^cosmos /^sdss* as a function of redshift". 



redshift 


type 


P< 


P> 


0.3 


type 2 


1.386 ± 0.327 


1.664 ± 0.513 


0.5 


type 2 


0.710 ± 0.163 


0.686 ± 0.214 


0.7 


type 2 


0.960 ± 0.196 


0.681 ± 0.197 


0.9 


type 2 


0.920 ± 0.195 


0.616 ± 0.185 


0.3 


type 2.3 


0.904 ± 0.398 


0.798 ± 0.502 


0.5 


type 2.3 


0.710 ± 0.274 


0.610 ± 0.336 


0.7 


type 2.3 


1.169 ± 0.393 


0.706 ± 0.349 


0.9 


type 2.3 


1.550 ± 0.528 


1.142 ± 0.551 


0.3 


type 2.2 


1.639 ± 0.585 


2.379 ± 1.241 


0.5 


type 2.2 


0.737 ± 0.262 


0.819 ± 0.445 


0.7 


type 2.2 


0.832 ± 0.270 


0.811 ± 0.411 


0.9 


type 2.2 


0.632 ± 0.215 


0.372 ± 0.205 


0.3 


type 2.1 


1.954 ± 1.392 


1.311 ± 0.934 


0.5 


type 2.1 


0.792 ± 0.573 


0.582 ± 0.415 


0.7 


type 2.1 


1.089 ± 0.706 


0.446 ± 0.305 


0.9 


type 2.1 


0.707 ± 0.473 


0.311 ± 0.217 


0.3 


type 2.0 


1.473 ± 0.903 


2.651 ± 1.973 


0.5 


type 2.0 


0.578 ± 0.364 


0.724 ± 0.584 


0.7 


type 2.0 


0.725 ± 0.403 


0.681 ± 0.513 


0.9 


type 2.0 


0.338 ± 0.209 


0.345 ± 0.279 


Relative densities 


P = ^COSMOS 


/^SDSS* for disk 



galaxies as a function of redshift and S/D-category. Here 
^COSMOS ^SDSS* denominate the comoving num- 
ber density per unit size as reported in the plots of the size 
function in Figures 13 and 14 for the COSMOS sample 

in the redshift range z G [0.6, 0.8[, and in Figures 15 and 
17 for the artificially redshiftcd and brightened sample of 
the SDSS galaxies. p< stands for the relative density of 
disks with intermediate sizes of ri/2 G [5 kpc,7 kpc[ and 
io> for that of large disks with ri/2 G P kpc, 17 kpc[. 
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Fig. 1. — Selection of objects with high values of Xgim2D but satisfactory fits to the smooth 
component of the surface brightness profile of the galaxy. Prom left to right, the figure 
shows the original stamp image, the cleaned image, the GIM2D model and the residual 
image (c.f. section 3). For each object, the bar in the original galaxy image indicates a scale 
of 2 arcseconds. An intensity scale for the normalized residual images is reported on their 
right-hand side. Model parameters are hsted in Table 1. 
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Fig. 2. — Original images, cleaned images, GIM2D model and residual image for three 
galaxies of ZEST morphological type 1. Color and size scales are as in Figure 1; the values 
of the physical parameters are listed in Table 1. 
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Fig. 3. — Original images, cleaned images, GIM2D model and residual image for three 
galaxies of ZEST morphological type 2.3. Color and size scales are as in Figure 1; the values 
of the physical parameters are listed in Table 1. 



Fig. 4. — Original images, cleaned images, GIM2D model and residual image for three 
galaxies of ZEST morphological type 2.2. Color and size scales are as in Figure 1; the values 
of the physical parameters are listed in Table 1. 
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Fig. 5. — Original images, cleaned images, GIM2D model and residual image for three 
galaxies of ZEST morphological type 2.1. Color and size scales are as in Figure 1; the values 
of the physical parameters are listed in Table 1. 
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Fig. 6. — Original images, cleaned images, GIM2D model and residual image for three 
galaxies of ZEST morphological type 2.0. Color and size scales are as in Figure 1; the values 
of the physical parameters are listed in Table 1. 
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Fig. 7. — Original images, cleaned images, GIM2D model and residual image for three 
galaxies of ZEST morphological type 3. Color and size scales are as in Figure 1; the values 
of the physical parameters are listed in Table 1. Due to the irregularity of the distribution 
of light in this class of objects, these GIM2D fits are inevitably the most unrehable. 
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n n 

Fig. 8. — The comparison of the non-parametric diagnostics A (asymmetry), C (concentra- 
tion), M20 (second order moment of the brightest 20% of the pixels) and G (Gini coefficient) 
with the Sersic index n. The regular trends and tight relations in these plots support the 
reliability of the GIM2D fits. The clustering of points at n = 9 refiects our restriction of the 
Sersic index to the interval n e [0.2, 9]. 
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Fig. 9. — Left half: Distribution of Scrsic indices for galaxies classified by ZEST as type 1 
(early-types), type 2 (disk galaxies) and type 3 (irregular and peculiar galaxies). Right half: 
Distribution of Sersic indices for sources ranging from purely exponential disks (type 2.3) to 
bulge-dominated disks (type 2.0). 
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""i/z [kpc] r,/2 [kpc] 

Fig. 10. — Left half: Distribution of half-light radii (measured in kpc) for galaxies classified 
by ZEST as type 1, type 2 and type 3. Right half: Distribution of half-light radii for objects 
ranging between pure exponential disks (type 2.3) and strongly bulge-dominated objects 
(type 2.0). 
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redshift 



Fig. 11.— The distribution of all Cycle 12 COSMOS galaxies with / e [16, 22.5] and absolute 
S-band magnitude Mb at redshifts z < 1.5. 




1 10 1 10 1 10 1 10 

Fig. 12. — Distribution of absolute S-band magnitude Mb vs. physical size for type 2 disk 
galaxies as a function of redsliift. Horizontal lines show the average magnitude cuts applied 
in the derivation of the size function within the different redshift bins (c.f. Section 4.2) and 
inclined dotted lines are lines of constant surface brightness. 
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Fig. 13. — Evolution of the size function between redshift 0.2 and 1 for the entire sample 
of disk galaxies. Four different bins of redshift, centred at 2; ~ 0.3, 0.5, 0.7 and 0.9, are 
displayed. The vertical error bars reproduce the Poisson errors. Areas shaded in grey 
indicate the errorbar on the size function induced by random errors in redshift, size and 
observed magnitude (c.f. Section 4.3). The solid and dashed lines in black show fits to the 
relation log{^) — a • ri/2 + P above 5 kpc, calculated for two sets of ZEBRA photometric 
redshifts obtained without and with corrections to the photometric catalogs, respectively 
(c.f. Section 4.1; Table 4 lists the best fit parameters of the solid lines). 
Apart from a modest steepening of the slope with increasing redshift, the size function of 
disk galaxies is roughly constant in the range 0.2 < z < 1.0. 
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Fig. 14. — The evolution of the size function shown individually for all sub-classes of the 
disk galaxy population in COSMOS. The sub-classes are displayed by increasing S/D-ratio 
from top to bottom. The errors and the solid and dashed black hues are as in Figure 13. 
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Fig. 15. — Normalized distribution of absolute B-band magnitudes (top), Sersic indices 
(centre) and half-light radii (bottom) in the sample of artificially redshifted SDSS galax- 
ies (black histogram; brightened sample - grey histogram; unbrightened sample) and the 
COSMOS data set (red histograms) in the range z G [0.6,0.8[. 
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Fig. 16. — Size function for COSMOS disk galaxies as plotted in Figure 13 (black squares) in 
the range z G [0.6, 0.8 [ and for brightened (filled grey points) and unbrightened (empty grey 
points) artificially redshifted SDSS galaxies. The data points of the SDSS size functions have 
been slightly offset from their nominal value of ri/2 for better legibility. Systematic effects 
on the GIM2D size measurements are identical to those discussed for the COSMOS sample 
(see caption Figure 13, and Appendix Al). 

The COSMOS size function in the range z e [0.6, 0.8[ and that of the simulated SDSS* 
galaxies closely resemble each other. The steeper slope of the size function in COSMOS 

relative to that of artificially redshifted SDSS galaxies confirms the tendency observed in 
Figures 13 and 14, implying a small deficit of the largest disks (ri/2 > 7 kpc) in the COSMOS 
sample compared to the local universe. 
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Fig. 17. — Size function for COSMOS disk galaxies (reported with black squares and split by 
relative importance of the bulge component) in the range z G [0.6, 0.8 [ (as plotted in Figure 
14) and for brightened (filled grey points) and unbrightened (empty grey points) artificially 
redshiftcd SDSS galaxies. 

Broadly speaking, a general constancy of the size function over half the Hubble time equally 
is seen for galaxies of all bulge-to-disk ratios. However, the four panels above indicate 
that the steepening of the slope of the size function is more pronounced for galaxies with 
measurable bulges than for bulgeless disks. Therefore, the deficit of large disks at earher 
epochs is apparently caused by a deficit of bulged disk galaxies at high redshift. 
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Fig. 18. — The density of COSMOS disk galaxies with half-hght radii in the range 5kpc- 
7kpc and >7kpc, normahzed to the corresponding brightened SDSS sample, as a function of 
redshift (c.f. Table 6 for the values of the relative density p — ^cosmos /^sdss*) ■ The point 
at ^ = (equal to 1 by definition since p = ^sdss*/^sdss* ) is plotted to show the error on 
the z = measurement based on the SDSS sample. The evolution of the SDSS-normalized 
density is not affected by the systematic errors in the GIM2D estimates of the half-light 
radii discussed in Appendix Al. Black and grey symbols refer to the two sets of photometric 
redshifts from ZEBRA, obtained without and with correction to the photometric catalogs 
(c.f. Section 4.1 and Figures 13 and 14). Uncertainties in the estimates of the photometric 
redshifts do not change our conclusions. 
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Fig. 19. — As in Figure 18, but shown individually for disk galaxies with different B/D-reitio. 
Black and grey symbols too are as in Figure 18. 

A constant (or possibly even slightly higher) number density of large bulgeless disks is ob- 
served at ^ ~ 1, while disk galaxies with a prominent bulge are three time less frequent than 
in the local universe at that redshift. 
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Fig. 20. — Dependence of the error in the recovered Sersic index n on input magnitude 
Iabau- On the left we plot the distribution of normalized difference between the nominal n 
value and the Sersic index recovered by GIM2D for nj„ = 1 galaxies, on the right the same 
for objects with a nominal Sersic index of n^^ = 4. 
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Fig. 21. — Upper panel: Distribution of normalized difference between the nominal input 
half-light radius of the simulated galaxies and the half-light radii recovered with the GIM2D 
fits as a function of magnitude. Lower panel: As above, but showing the distribution for 
galaxies with input half-light radii smaller and larger than 10 pixels. 
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Fig. 22. — Variation of GIM2D output for a specific galaxy when fit using 250 different point 
spread functions from the collection of PSFs employed in our work. The PSFs have been 
derived for different individual ACS tiles, and vary as a function of position within a single 
tile. 
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